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Abstract 

■^jp '. The mixing of the known quarks with new heavy singlets can modify signif- 



o 
o 



icantly some observables with respect to the Standard Model predictions. We 



CN , analyse the range of deviations permitted by the constraints from precision elec- 

troweak data and flavour-changing neutral processes at low energies. We study 
top charged current and neutral current couplings, which will be directly tested 
at top factories, like LHC and TESLA. We discuss some examples of observables 
from K and B physics, as the branching ratio of Kl — ► ir®vu, the B® mass dif- 
ference or the time-dependent CP asymmetry in B® — ► DfD~, which can also 



show large new effects. 



(N : 1 Introduction 
o 

The successful operation of LEP and SLD in the past few years has provided pre- 
Qh| cise experimental data [TJ E] with which the Standard Model (SM) and its proposed 

extensions must be confronted. The results for e'/e have converged providing 
evidence for direct CP violation in the neutral kaon system. In addition, B factories 
^ ■ have begun producing data on B decays and CP violation, which test the Cabibbo- 

Kobayashi-Maskawa (CKM) matrix elements involving the top quark and the CKM 
phase 5. However, the determination of most parameters involving the top quark is still 
strongly model-dependent. While the CKM matrix elements that mix light quarks are 
extracted from tree-level processes (and hence their measurement is model-independent 
to a large extent), the charged current couplings V t d and V ts are derived from one-loop 
processes 0, to which new physics beyond the SM may well contribute. The Teva- 
tron determination of V t b in top pair production jB] is obtained assuming 3x3 CKM 
unitarity, and the neutral current interactions of the top with the Z boson remain vir- 
tually unknown from the experimental point of view. This fact contrasts with the high 
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precision achieved for the couplings of the b and c quarks at LEP and SLD, obtained 
from the ratios Rb, R c and the forward-backward (FB) asymmetries A F #, A F #. 

The situation concerning CP violating phases is better (see for instance Ref. |Z1E])- 
Few years ago, the single phase S present in the CKM matrix could merely be adjusted 
to reproduce the experimental value of the only CP violation observable available, e 
in the kaon system. With the resolution of the conflict between the NA31 and E731 
values of e'/e, and the recent measurement of the CP asymmetry a^K s m the B system 
QDj, there are two new CP violation observables, both in agreement with the SM 
predictions, which allow to test the CKM picture of CP violation. Experiments under 
way at B factories keep investigating other CP asymmetries to dig out the phase 
structure of the CKM matrix. Likewise, the knowledge of the top quark properties will 
improve in the next years, with the arrival of top factories, LHC and TESLA |TT) I12j. 
For instance, single top production at LHC [T3J EH UH] will yield a measurement of V t b 
with an accuracy of ± 7%. In top pair production, the angular distributions of the top 
decay products will provide a very precise determination of the structure of the Wtb 
vertex, even sensitive to QCD corrections [TKj . The prospects for V t d and V ts are less 
optimistic due to the difficulty in tagging the light quark jets. 

Before top factories come into operation, it is natural to ask ourselves how large 
the departures from the SM predictions might be. Answering this question means 
knowing how precisely one can indirectly fix the allowed values of the least known 
parameters, taking into account all the present relevant data from electroweak precision 
measurements and from kaon, D and B physics. We will show that there is still 
large room for new physics, which may manifest itself in the form of deviations of the 
properties of the known quarks from SM expectations. This is especially the case for 
the top quark, whose couplings are poorly known, and also for rare K decays and CP 
asymmetries in the B systems, which are currently being probed at K and B factories. 

With this aim we study a class of SM extensions in which Q = 2/3 up- type or Q = 
-1/3 down-type quark singlets are added to the three SM families [T7 j lT8 | rm | 120 1 121 ] l22 | 
EH 123 EH1 EH EZ| • These exotic quarks, often called vector-like, have both their left and 
right components transforming as singlets under SU(2)^, and therefore their addition 
to the SM quark content does not spoil the cancellation of the triangle anomalies. In 
these models, which are described in the next Section, 3x3 CKM unitarity does not 
necessarily hold, and mixing of the new quarks with the standard ones can lead to 
sizeable departures from the SM predictions |2BJ EH EOl dJ] • For instance, the CKM 
matrix elements V t d, V t3 and V t b and the top neutral current couplings with the Z boson 
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can be quite different from SM expectations. The ratio of branching fractions of the 
"golden modes" By(Kl — > 7r ( Vz/)/Br(.ft' + — > 7i + uu) can have an enhancement of one 
order of magnitude with respect to the SM prediction, and the time-dependent CP 
asymmetry in the decay B® — ► D+D~, which is predicted to be very small in the SM, 
can vary between —0.4 and 0.4. 

Apart from their simplicity and the potentially large effects on experimental ob- 
servables, there are several theoretical reasons to consider quark isosinglets. Down 
singlets appear in grand unification theories [221 E21 EH] , for instance those based on 
the gauge group E 6 (in the 27 representation of E 6 a Q = —1/3 singlet is associated 
to each fermion family). The presence of down singlets does not spoil gauge coupling 
unification, as long as they are embedded within the 27 representation of Eq |23| I2fij. 
When added to the SM particle content, they can improve the convergence of the cou- 
plings, but not as well as in the minimal supersymmetric SM [33]. Models with large 
extra dimensions with for instance t R in the bulk predict the existence of a tower of 
Q = 2/3 singlets Tjj R - If there is multilocalisation the lightest one, T^/ R , can have a 
mass rriT ~ 300 GeV or larger and an observable mixing with the top [3S]- Similarly, 
if b R is in the bulk, there exists a tower of Bj^ R , but the mixing with SM fermions is 
suppressed in relation to the top case if the Higgs is restricted to live on the boundary. 

There are three recent studies regarding the constraints on models with extra sin- 
glets [Hni GUI EE] • In this work, we extend these analyses in three main aspects: (i) 
We include up singlets, as well as down singlets, referring to them as Model I and E, 
respectively, (ii) We study the limits on V t d, V ts , V t b, top neutral current couplings and 
other observables not previously analysed; (m) We take a larger set of experimental 
constraints into account: the correlated measurement of R^, R c , Ap B , Apg, Ab, A c ; 
oblique corrections; the |5j7ig|, l^mej arid |5m£>| mass differences; the CP violation 
observables e, e'/e, a^K s ] the decays b — > 57, b — > s[i + fi~, b — > se + e~, K + — > -n + vu 
and Kl — > /i + /i~; vN processes and atomic parity violation. In addition, we examine 
several other potential restrictions, which turn out to be less important than the pre- 
vious ones. We allow mixing of all the generations with either Q = 2/3 or Q = —1/3 
exotic quarks, and we consider that one or two singlets can mix significantly, though 
for brevity in the notation we always refer to one extra singlet T or B. 

This paper is organised as follows: In Section [2] the main features of the models are 
described. In Section El we summarise the direct limits on CKM matrix elements and 
the masses of the new quarks. In Section 0] we review the constraints from precision 
electroweak data: R b , R c , asymmetries and oblique corrections. In Section E] we focus 
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our attention on flavour- changing neutral (FCN) processes at low energies: meson 
mixing, B decays and kaon decays. The various constraints on the Z couplings of the 
u, d quarks are studied in Section El We introduce the formalism necessary for the 
discussion of some observables from K and B physics in Section We present the 
results in SectionlHl and in Section|H]we draw our conclusions. In Appendix|X]we collect 
the common input parameters for our calculations, and in Appendix El the Inami-Lim 
functions needed. The statistical prescriptions used in our analysis are explained in 
Appendix O 

2 Brief description of the models 

In order to fix our notation briefly, in this Section we will be a little more general than 
needed in the rest of the paper (see for instance Ref. f^lj for an extended discussion 
including isodoublets and mirror quarks too). We consider a SM extension with N 
standard quark families and n u up, rid down vector-like singlets. The total numbers of 
up and down quarks, M u = N + n u and Md = N + rid, respectively, are not necessarily 
equal. In these models, the charged and neutral current terms of the Lagrangian in 
the weak eigenstate basis can be written in matrix notation as 



with {u L , d L ) doublets under SU(2)l of dimension N in flavour space. These terms 
have the same structure as in the SM, with N generations of left-handed doublets in 
the isospin-related terms, but with all the M u , Md fields in Jem- The differences show 
up in the mass eigenstate basis, where the Lagrangian reads 



Here u = (u, c, t,T, . . .) and d = (d, s, b,B, . . .) are Af u and Md dimensional vectors, 
and X u , X d are matrices of dimension M u x M u , Md x Md, respectively. In general the 
M u x Md CKM matrix V is neither unitary nor square. 

The most distinctive feature of this class of models is the appearance of tree-level 
FCN couplings in the mass eigenstate basis, originated by the mixing of weak eigen- 
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states with the same chirality and different isospin. These FCN interactions mix left- 
handed fields, and are determined by the off-diagonal entries in the matrices X u , X d . 
On the other hand, the diagonal Zqq terms of up- or down-type mass eigenstates q are 
(dropping here the superscript on the X matrices) 

C L = ^Xgg — 2Q q S w , 

c q R = —2Qqsly , (3) 

with the plus (minus) sign for up (down) quarks. With these definitions, the flavour- 
diagonal Zqq vertices read 

Czqq = -^^^L + qRl^qn) . (4) 

For a SM-like mass eigenstate without any left-handed singlet component, X qq = 1, 
X qq i = for q' ^ q, and it has standard interactions with the Z boson. For a mass 
eigenstate with singlet components, < X qq < 1, what implies nonzero FCN couplings 
Xq q i as well. 

Let us write the unitary transformations between the mass and weak interaction 
eigenstates, 



u° L = U uL u L , u° R = U uR u R , 

d° L = U dL d L , d° R = U dR d R , (5) 



where U qL and U qR are M q x M q unitary matrices. The weak interaction eigenstates 
q\ R include doublets and singlets. It follows from Eqs. (jlf2|) that 



/uL\* i idL 



ru (1 juL\* 7 Jl 



Kb = {u^rK, x d aT = (uf a L yuf T L (6) 



with i, j running over the left-handed doublets, a, j3 = u,c,t,T, . . . and <t,t = d, s, b,B, . . 
From these equations it is straightforward to obtain the relations 

x u = vv\ 

X d = V^V, (7) 

and to observe that X u = (X u y, X d = (X d y. Furthermore, we can see that in general 
V is not an unitary matrix. We will restrict our discussion to models where either 
rid = or n u = 0, i. e. we will consider either up singlets or down singlets, but not 



5 



both at the same time. In this context V is a submatrix of a unitary matrix, and in 
each case we can write 

i=l i=N+l 
N Af d 

xj T = y £Kv*r = s„- E Kv*. (8) 

i=l i=N+l 

It is enlightening to observe that for a ^ j3, a ^ r, we have the inequalities [351 I4*U] 

\X d ar \ 2 < (1 - Xt){l - X d TT ) . (9) 

Therefore, if for instance X£ a = 1 (that is, if the diagonal Z vertex is the same as in 
the SM) the off-diagonal couplings involving the quark a vanish. As a rule of thumb, 
FCN couplings arise at the expense of decreasing the diagonal ones. This fact has 
strong implications on the limits on FCN couplings, as we will later see. 

The equality for X u in Eq. © holds in particular if n u — 1. Likewise, the equality 
for X d holds when rid = 1- This implies that the introduction of only one extra 
singlet mixing significantly (as it is usually done in the literature) results in additional 
restrictions in the parameter space, and in principle may lead to different predictions. 
Moreover, for n u = 1 or = 1 the CKM matrix has three independent CP violating 
phases, whereas for n u = 2 or = 2 there are five. Hence, in our numerical analysis 
we will consider also the situation when two singlets can have large mixing, n u = 2, 
nd = or n u = 0, nd = 2, to give a more complete picture. In the rest of the paper we 
write the expressions for only one extra singlet for simplicity. 



3 Direct limits 

Even though in these SM extensions the 3x3 CKM matrix is not unitary, in the two 
examples under study it is still a submatrix of a 4 x 4 unitary matrix V. The direct 
determination of the moduli of CKM matrix elements [S] in Tabled not only sets direct 
limits on these CKM elements themselves but also unitarity bounds on the rest. After 
the requirement of V t b ~ 1 from precision electroweak data (see Section HJ) these bounds 
are stronger. In this case, the Tevatron constraint [Oj 

0.97i°Ji (10) 





v tb \ 


2 


\Vu\ 


2 + \ 




| 2 + 


v tb \ 


2 



6 



v ud \ 


0.9735 ±0.0008 


v us \ 


0.2196 ±0.0023 


v ub \ 


(3.6 ± 1.0) x 10- 


Veil 


0.224 ±0.016 


v cs \ 


0.97 ±0.11 


v cb \ 


0.0402 ±0.0019 



Table 1: Direct measurements of CKM matrix elements. V ub is obtained from \V cb \ and 
the ratio \V ub /V cb \. 

is automatically satisfied. 

The non-observation of top decays t — > qZ, q = u, c at Tevatron ^T] provided the 
first limit on FCN couplings involving the top, \X qt \ < 0.84 (from now on we omit 
the superscript when it is obvious). These figures have improved with the analysis of 
single top production at LEP in the process e + e~ — > tq + iq, which sets the bounds 
\X qt \ < 0.41 [42J. LEP limits are model-dependent because in single top production 
there might possibly be contributions from a jtq effective coupling. These vertices are 
very small in most SM extensions, in particular in models with quark singlets [13] , thus 
in our case the photon contribution may be safely ignored. 

As long as new quarks have not been observed at Tevatron nor LEP, there are 
various direct limits on their masses, depending on the decay channel analysed [5] . We 
assume m T , m B > 200 GeV in our evaluations. 

4 Limits from precision electroweak data 

4.1 Rb, R c and FB asymmetries 

In the discussion after Eqs. (jOJ) we have observed that FCN interactions can be bounded 
by examining the deviation from unity of the diagonal ones. This is a particular 
example of a more general feature of these models, that the isosinglet component of 
a mass eigenstate can be determined from its diagonal couplings with the Z boson. 
In this Section we will explain how the experimental knowledge of R b , R c and the FB 
asymmetries of the b and c quarks constrains their mixing with isosinglets. We will 
study in detail the case of the bottom quark; the discussion for the charm is rather 
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alike. 

Rb is defined as the ratio 

Ro _ r(z -> bl) 

r(Z — > hadrons) 

The partial width to hadrons includes uu, dd, ss, cc and bb. The numerator of this 
expression is proportional to {c b L ) 2 + {c b R ) 2 plus a smaller term proportional to mj. The 
pole FB asymmetry of the b quark is defined as 

,0,6 tx(cos^ > 0) - a(cosfl < 0) 

a(cos > 0) + ^(cos < 0) ' 1 j 

where 9 is the angle between the bottom and the electron momenta in the centre of 
mass frame . The coupling parameter Ab of the bottom is defined as 

A _ ( C Z,) 2 ~ ( C R) 2 /lO\ 

b ~ (4) 2 + & ■ {16} 

It is obtained from the left-right-forward-backward asymmetry of the b quark at SLD, 
and considered as an independent parameter in the fits, despite the fact that the FB 
asymmetry can be expressed as A F ^ = 3/AA e Ab, with A e the coupling parameter of 
the electron. 

At tree-level, c b L = —Xbb + 2/3s 2 y , c 9 R = 2/3s 2 y, hence in a first approximation the 
mixing of the b quark with down singlets in Model E decreases Xbb from unity and thus 
decreases Rb, Ab and A F #. The effect of some electroweak radiative corrections can be 
taken into account using an MS definition of the sine of the weak angle, s\ = 0.23113 
[3] and for the electron coupling an "effective" leptonic sin 2 #^ t = 0.23137 pQ. Other 
electroweak and QCD corrections that cannot be absorbed into these definitions are 
included as well They are of order 0.6% for u, c and —0.25 % for d, s, 

b. Furthermore, for the bottom quark there is an important correction originated by 
triangle diagrams involving the top 



5c b L = 2^\V tb \ 2 F(x t ) (14) 

(note that we use a different normalisation with respect to Ref. |46J, with x t = 
(m t /M w ) 2 and 



F(x t ) = —y- [x t + 2.880 log x t - 6. 716 + (8.368 log x t -3A08)/x t 

+ (9.126 log x t + 2.260)/x 2 + (4.043 log x t + 7.410)/x t 3 + . . .] . (15) 



1 These two observables do not include the photon contributions, and A FB is defined for mass- 
less external particles. They are extracted from the experimental measurement of e + e~ — > bb after 
correcting for the photon contribution, external masses and other effects G] ■ 
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We have omitted the imaginary part of F(x t ) since it does not contribute to 5c b L . This 
large correction ~ {m t /M w ) 2 is a consequence of the non-decoupling behaviour of the 
top quark, and CKM suppression makes it relevant only for the bottom. It decreases 
the value of Rb by 4cr and has the indirect effect of increasing R c slightly. Its inclusion 
is then crucial to compare the theoretical calculation with experiment. In Table El 
we collect our SM predictions for R b , R c , A FB , Apg, Ab and A c calculated using the 
parameters in Appendix El together with the experimental values found in Ref. pQ. 
The masses used are MS masses at the scale Mz- The correlation matrix necessary for 
the fit is in Table El 





SM 


Experimental 


Total 




prediction 


measurement 


error 


Rb 


0.21558 


0.21646 


0.00065 


R c 


0.1722 


0.1719 


0.0031 


A o,b 

^FB 


0.1039 


0.0990 


0.0017 


A 0,c 
^FB 


0.0744 


0.0685 


0.0034 


A b 


0.935 


0.922 


0.020 




0.669 


0.670 


0.026 



Table 2: SM calculation of Rb, R c , A F ^, A FB , Ab, A c and experimental values. 

Rb Rc ^FB ^FB ^b A c 

R h 1.00 -0.14 -0.08 0.01 -0.08 0.04 

R c -0.14 1.00 0.04 -0.01 0.03 -0.05 

A°'b -0.08 0.04 1.00 0.15 0.02 0.00 

A% 0.01 -0.01 0.15 1.00 0.00 0.01 

A b -0.08 0.03 0.02 0.00 1.00 0.13 

A c 0.04 -0.05 0.00 0.01 0.13 1.00 



Table 3: Correlation matrix for the experimental measurements of Rb, R c , A F '^, A FB , 
Ab and A c . 

The mixing of the b quark with down isosinglets decreases Vtb, making this negative 
correction smaller in modulus. This is however less important than the effect of the 
deviation of X&& from unity. The net effect is that in Model II X^, and hence also V t b, 
are tightly constrained by R b to be very close to unity. 
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In Model I the mixing of the top with singlets modifies the Ztt interactions, and 
the expression for Sc b L in Eq. (JT3j) must be corrected accordingly (see Ref. [H] and 
also Ref. [2E|). The decrease in X u can be taken into account with the substitution 
F -»• F + F 2 , with 

F 2 (x t ) = ( 2 - ■ ( 16 ) 

Moreover, there are additional triangle diagrams with the mass eigenstate T replacing 
the top, or involving t and T. The T quark contribution is added to Eq. (jl4j) as the 
top term but multiplied by | Vrb\ 2 - The t — T contribution is given by V t * b Vrb F 3 (x t , xt), 
with x T = (m T /M w ) 2 and 2 

1 ReX tT 



F 3 (x t ,x T ) 



2s 2 



Xt 2 xf N 
■ log xt log x t 



xt — x t \xt — 1 x t — 1 



x t x T ( x T x t 

log Xt r log x t 



xt — x t \xt — 1 £ t — 1 



(17) 



In Model I this radiative correction gives the leading effect on Rb of the mixing. 
However, the presence of the new quark may make up for the difference in the top 
contribution. Should the new mass eigenstate be degenerate with the top, uit = m t 
and Xt = x t , one can verify that 

\V tb \ 2 F(x t ) + \V Tb \ 2 F(x T ) = \V tb \ 2 [F(x t ) + F 2 (x t )] + \V Tb \ 2 [F(x T ) + F 2 {x T )\ 

+V t tV Tb F 3 (x t ,x T ), (18) 

as intuitively might be expected. Since (|Vt;,| 2 + |Vt6| 2 ) = |Vt6|§M; ^ ms means that for 
degenerate t, T the correction has the same value as in the SM (and in this situation the 
terms with F 2 and F 3 cancel each other). For itlt ~ m t , 5c b L has a similar magnitude as 
in the SM and low values Vtb ~ 0.6 are allowed. For heavier T, the size of this radiative 
correction sets limits on the CKM angle Vrb, and thus on Vtb- 

The study of the charm mixing and the constraints on its couplings from R c , A c 
and Apg is completely analogous (interchanging the role of up and down singlets). In 
principle, the presence of a new heavy down quark B induces a large m^-dependent 
correction, but this is suppressed by the CKM factor |V^b| 2 and hence the analysis is 
simplified. The pole FB asymmetry of the quark s has also been measured recently, 



2 In obtaining Eq. i|17|) from the results quoted in Ref. we have assumed a CKM parameterisation 
with V t * b Vrb real. This is our case with the parameterisations used in the numerical analysis in 
Section 0D 



10 



Apg = 0.1008 ± 0.0120 j3H], though not nearly with the same precision as the b and c 
asymmetries. This determination assumes that the FB asymmetries of the u, d quarks 
and the Z branching ratios are fixed at their SM values and thus cannot be properly 
taken as a direct measurement. We do not include it as a constraint, but anyway 
we have checked that at this level of experimental precision it would not provide any 
additional constraint on the model. 



4.2 Oblique parameters 

The oblique parameters S, T and U [SHI EI| are used to summarise the effects of 
new particles in weak currents in a compact form. Provided these particles are heavy 
and couple weakly to the known fermions, their leading effects in processes with only 
SM external particles are radiative corrections given by vacuum polarization diagrams 
(oblique corrections), rather than triangle and box diagrams (direct corrections). We 
will use the definitions ESj 

n 3 y(Ml)-n 3 y(o) 

S = ~ 16n Ml ' 

t = 77^^1(0)^33(0)], 

Ivl Z b W c W 

/ - n 11 (A4.)-n 11 (o) iUMp-iuo^ 

U = 161 { Ml Ml j ■ (19) 

They are equivalent to the ones used in Ref. (Hj, as can be seen by a change of basis. In 
these expressions only the contributions of new particles are meant to be included. Ra- 
diative corrections from SM particles must be treated separately because their leading 
effects are direct, not oblique. The parameters S, T, U are extracted from precision 
electroweak observables, and their most recent values are in Table |U The contributions 
to S, T and U of an arbitrary number of families plus vector-like singlets and doublets 
have been computed in Ref. [231 • m our models there are no exotic vector-like doublets, 
hence right-handed currents are absent and their expressions simplify to 



I M u Af„ M d 

\q=1ct=1 (3<a T<er 



/ N u M d A/-„ Md 



N I Md AT„ M d 

u = EEi^i 2 ^+(^'^)-Ei^i 2 ^(y-^)-Ei x -i 2 x+(y-^)] , (20) 

\a=lcr=l f3<a T<a 
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where N c = 3 is the number of colours, y,i = (rrii/Mz) 2 and we use the MS definition 
of s^y, as well as MS masses at the scale M Z - The functions multiplying the mixing 
angles are 



V± (2/1,2/2) 

#+(2/1,2/2) 
x+ (2/1, 2/2) 



22 Vl + Uy 2 
9 



22/12/2 , 2/1 

2/i ± 2/2 log — , 

2/i - 2/2 2/2 

yi + V2 _ (yi - V2) 2 



g log — + — — — /(yi, 2/1) + — — /(S/2, 2/2) , 



2 

J/1-1 



(l/i - 2/2 ) : " 
6 



7(2/1,2/1) + ^-z—f(y2,y2) 



y\ + 2/1 
2(2/1 - 2/2) 



log 



2/1 
2/2 



+ 



6 - / 6 

1 _ 2/1 +2/2 _ (2/1 ~ 2/2) 2 
3 



7(2/1,2/2) • 



6 6 
The function / is defined as 

-2y/A (arctan yi ^ +1 - arctan 

^ZAlog «i+«2-i+^ 
S/1+1/2-1-V-A 



(21) 



/(2/i,2/2) 



A > 



A < 



(22) 



with A = — 1 — — y\ + 2j/i + 2y 2 + 2yij/2- The functions ■?/>, 0, x are symmetric under 
the interchange of their variables, and 9, x satisfy 9(y,y) = 0, xiViV) — 0. 



S -0.03 ±0.11 
T -0.02 ±0.13 
U 0.24 ±0.13 



Table 4: Experimental values of the oblique parameters. 



These expressions are far from transparent, and to have a better understanding of 
them we will examine the example of an up singlet mixing exclusively with the top. In 
this limit, the new contributions are 

S = ^{\V Tb \ 2 [ip + (y T ,y b ) -ip + (yt,yb)]-\X T t\ 2 ip+(yT,yt)} , 

T = 1fi 2 2 {\V Tb \ 2 [9 + (yT,y b )-9 + (y t ,y b )]-\X Tt \ 2 9 + (y T ,y t )} , 

u = -^{\v Tb \ 2 lx+(yT,y b ) -x+(yt,y b )]-\x T t\ 2 x+(yT,yt)} ■ (23) 
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The factors |Vr&| 2 , |^t*| 2 describing the mixing of the quark T are not independent: 
as can be seen from the results in Section EJ for n u = 1 we have the relation 



IX 



Tt\ 



\v Tb v tb \ 2 



IV, 



Tb\ 



;i-i^i 2 ) 



(24) 



For t and T degenerate, T and U would automatically vanish independently of the 
mixing, and S = —0.16 |Xj^| 2 . In order to obtain a simple approximate formula when 
m T 3> m t , we approximate 1 — \Vrb\ 2 ~ 1 and keep only the leading order in y T . 



(Needless to say, we use Eqs. (J20H22J) for our fits. 
yt, this yields 



Using the numerical values of y b , 



S 



T 



U 



Nc \v I 

Nc \V ' 2 



2tt' 



'Tb\ 



-0.34 + 0{vt 1 ) 

\V Tb \ 2 [-18.4 + 7.8 log y T + 0{y T l ) 
0.60 + 0(y T l 



(25) 



These expressions give a fair estimate of the effect of the top mixing in the oblique 
parameters. We notice that the effect on S, U is very small, S = —0.16 | V^rf, | 2 , U = 
0.29 \V Tb \ 2 , but sizeable for T (for instance, T = 2.7 \V Tb \ 2 for m T = 500 GeV). Indeed, 
the T parameter bounds the CKM matrix element Vr b (and hence Vt b ) more effectively 
than the radiative correction to Rf, and better than low energy observables. 

In Model I the analysis is similar, but the constraints from Rb and FCN processes 
at low energies are much more restrictive than these from oblique corrections. 



5 Limits from FCN processes at low energies 

In this Section we discuss low energy processes involving meson mixing and/or decays. 
An important point is that almost all the observables analysed receive short-distance 
contributions from box and/or penguin diagrams with Q = 2/3 quark loops (otherwise 
it will be indicated explicitly). The top amplitudes are specially relevant due to the 
large top mass, and are proportional to V t * d V ts , V t * d V tb or V t *V tb (or their squares), de- 
pending on the meson considered. The observables are then sensitive to Vtd and Vt s . 
(Also to Vt b , but the most important restrictions on its modulus come from precision 
measurements examined in last Section.) Additionally, there are extra contributions 
in the models under study: either new box and penguin diagrams with an internal T 
quark in Model I or diagrams with tree-level flavour-changing neutral currents (FCNC) 
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mediated by the Z boson in Model E. In any case, the new terms depend on products 
of two elements of the fourth row of V (V£ d V T8 , V^ d V Tb or V^ s V Tb in Model I and FCN 
couplings Xds, Xdb or X s f> in Model I). 

The a priori unknown top and new physics terms are added coherently in the 
expressions of all these observables. Then, in principle there may exist a "conspiracy" 
between top and new physics contributions, with the first very different from the SM 
prediction and new physics making up for the difference. As long as we use a sufficiently 
exhaustive set of low energy observables and reproduce their experimental values, this 
possibility can be limited. This is because the products V td V ts , . . . , V^d^Ts or X ds , 
etc. appear in the expressions of these observables in combinations with different 
coefficients. 

Our observables for Models I and II include \5rriB s \, £, £'/£, the branching 

ratios for b — > se + e~, b — > s[i + fi~ , K + — > 7i + uu, Kl — > [i + and the CP asymmetry 
a^K s - For model I we use |<5m£>| as well. It must be stressed that they are all inde- 
pendent and give additional information that cannot be obtained from the rest. For 
example, if we remove e from the list we can find choices of parameters of our models 
for which all the remaining observables agree with experiment (the precise criteria of 
agreement used will be specified in Section |HJ) but e is more than 5 standard deviations 
from its measured value. This procedure applied to each one shows that none of them 
can be dismissed. 

Once the values of the observables in these sets are in agreement with experiment, 
the predictions for the mass difference and some other partial rates, like b —>■ 37, 

B — > suu, B — > B s — > /i + yU~, agree with SM expectations (b — > sj is nevertheless 

included in the fits). An important exception is the decay Kl — > 7r°z/z/, which will be 
studied in Section Several CP asymmetries can also differ from SM expectations, 
and are thus good places to search for departures from the SM or further restrict the 
models under consideration. 

In the rest of this Section we review the theoretical calculation within Models I and 
II of the observables listed above, together with their experimental status. 
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5.1 Neutral meson oscillations 



5.1.1 The AF = 2 effective Lagrangians 



The complete Lagrangian for Q = —1/3 external quarks in the presence of extra down 
singlets has been obtained in Ref. [54 3 and we follow their discussion except for small 
changes in the notation. We ignore QCD corrections for the moment and neglect 
external masses. The Lagrangians for K°, B° and B® oscillations are similar up to 
CKM factors, and for simplicity in the notation we refer to the kaon system. The box 
contributions can be written as 



/'box 



Gf a 



a/2 47ts 2 



a,(3=u,c,t 



(26) 



with = V* s V a d, etc. The function F is not gauge- invariant, and its expression in the 
't Hooft-Feynman gauge can be found e. g. in Ref. [S3]. (This and other Inami-Lim 
jSH] functions are collected in Appendix El) The terms involving the u quark can be 
eliminated using 



(27) 



and setting x u = 0, resulting in 
Gp a 



/'box 



y/2 47TS 2 



X! KdKd So( x a, x p) + 8X sd ^2 KdBo(x c 

a,f3=c,t 



Xsd 



a=c,t 



(28) 



where the gauge-independent function So is given in terms of the true box function F 
by 



S (x a ,x p ) = F{x a ,x f3 )-F{x a ,0)-F{0,x /3 )+F{0,0) 
and B Q is given in terms of F by 

4 5 (x Q ) = F(x Q ,0)-F(0,0). 



(29) 



(30) 



In addition there are two terms to be included in the Lagrangian. The first corresponds 
to Z tree-level FCNC, 



*-eff 



V2 



X 2 sd (s L Yd L ) (sl7/A) • 



(31) 
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The second is originated from diagrams with one tree-level FCN coupling and one 
triangle loop. Its contribution plus the B Q term in Eq. (|28j) can be compared with 
the short-distance effective Lagrangian for K° — ► (see Ref. jHl] for the details), 

concluding that the sum of both gives the gauge-invariant Inami-Lim function Yq with 
a minus sign. The full gauge-invariant AS = 2 effective Lagrangian then reads 



*-eff 



Gp 

71 



a 



O: 



X sd K 

a=c,t 



sd 



x isL-fdi) (s L ^d L ) . (32) 

We use the I superscript to refer to model E. The X% s term in Eq. (j2*6^1 is subleading 
with respect to C^ s , and it has been omitted. 

The SM Lagrangian can be readily recovered setting = in the above equation. 
In Model I the Lagrangian reduces to the SM-like box contributions but with terms 
involving the new mass eigenstate T, 



^eff — ~ 



Gp 
V2 



a 



E 



47FS W a,f3=c,t,T 



^sd^sd Soi X a, 



{sLl^dp) (s L 7aA) 



(33) 



In the B° and B° s systems the approximation of vanishing external masses is not jus- 
tified for the b quark. However, the two terms involving S (x c ) = S (x c ,x c ) and 
S (x c , x t ) are much smaller than the one with S (x t ) = S (x t , x t ) and can be neglected, 
and for the latter m& <C rrit and the approximation is valid. The effective Lagrangian 
for D° — D° mixing is more problematic and we will deal with it later. 

Short-distance QCD corrections are included in these Lagrangians as rj factors mul- 
tiplying each term in the usual way. These factors account for high energy QCD effects 
and renormalisation group (RG) evolution to lower scales When available, we use 
next-to-leading order (NLO) corrections [SH]. For the nonstandard contributions we 
use leading logarithmic (LL) RG evolution jSU]. The differences between LL and NLO 
corrections are minimal provided we use MS masses mj(mj) in the evaluations [H%j . 
Some representative QCD corrections for the new terms are 



Vtt 



[a s {m c )\ 
[oc s {m t )\ 



a s {m b ) 


6 

25 


'a s (M z )' 


23 


a s (m c )_ 




_ a s (m b ) _ 




a s (m T ) 


6 

21 







(34) 



a s (m t ) 

Here the superscripts K, B refer to the neutral mesons, and the subscripts to the term 
considered. 
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5.1.2 K° oscillations 

The element M 12 of the K° - K° mixing matrix is obtained from the effective La- 
grangian (see for instance Ref. |57|). 



K G 2 F M^f 2 K B K m K o 

12 ~ 127T 2 
c \ t K , 



(Ks) 2 V«S (x c ) + (\ ds f V %S (x t ) 



+2AS,A t ds ^5o(ar c ,z t )+Aj f J • (35) 

In this expression m^-o = 498 MeV is the K° mass, fx = 160 MeV the kaon decay 
constant taken from experiment and Bk = 0.86 ± 0.15 jOH| is the bag parameter. The 
QCD corrections are rj^ c = 1.38 ± 0.20, rj£ = 0.57, r/^ = 0.47 jHI] (we do not explicitly 
write the errors when they are negligible). The extra piece in models I and I is 

= (\1) 2 VttSo(x t ) + 2X c ds XlVc T S Q (x c , x T ) + 2\ t ds \ T dsV «S (x t , x T ) , 

^\ = -SX ds [X dsV §Y (x c ) + X^ t Y (x t )}+^^Xl, (36) 

with tjtt — 0.58, rj§ = 0.60. We estimate rj^ ~ rj^, ~ r/^ T , and expect that this is 
a good approximation because RG evolution is slower at larger scales. 

In the neutral kaon system the mass difference 5rriK can be written as 

5m K = 2 Re M& + 5m h £ (37) 

where the second term is a long-distance contribution that cannot be calculated reliably. 
For the first term we obtain (4.64 ± 0.68) x 10~ 3 ps" 1 within the SM, whereas = 
5.30 x 10~ 3 ps _1 . The large ~ 30 % long-distance contribution prevents us from using 
5rriK as a constraint on our models, but we observe anyway that the short- distance part 
2 Re Mf^ always takes values very close to the SM prediction once all other constraints 
are fulfilled. 

The CP violating parameter e is calculated as 3 

£ = e m/4 ^ — i 2 - (38) 
V2 5m K 

and in the SM it is close to its experimental value (2.282 ± 0.017) x 10" 3 after a 

proper choice of the CKM phase 5 (see Appendix EJ • The SM prediction with the 

3 This expression assumes a phase convention where V* d V us is real. For a rephasing-invariant 
definition of e see Ref. |62| . 
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phase 5 = 1.014 that best fits e, e'/e, a^ Ks and \5m B \ is e = (2.18 ± 0.38) x 10~ 3 . 
Notice that there is a large theoretical error in the calculation, mainly a consequence 
of the uncertainty in Bk, which results in a poor knowledge of the CKM phase that 
reproduces e within the SM. In Models I and I this parameter receives contributions 
from several CP violating phases and thus it cannot be used to extract one in particular. 
Instead, we let the phases arbitrary and require that the prediction for e agrees with 
experiment. 

5.1.3 B° oscillations 

The element M? 2 of the B° - B° mixing; matrix is 

with m B o = 5.279 GeV. We use f B = 200 ± 30 MeV, B B = 1.30 ± 0.18 from lattice 
calculations |B3]. The terms corresponding to Sq(x c ) and Sq(x c , x t ) have been discarded 
as usual, because in the SM, as well as in our models, they are numerically 2 — 3 orders 
of magnitude smaller than the So(x t ) term (the CKM angles are of the same order 
and the So functions are much smaller). The QCD correction is r)f t = 0.55. The 
nonstandard contributions are 

= (Xl) 2 V?TSo(x T ) + 2X t db X T db r ] f T S (x t ,x T ) } 

A| = -%X db \ dbV ?Mx t ) + ^ V B z Xl. (40) 

a 

The terms Sq(x c , xt) in A^ and Y (x c ) in A^ have been dropped with the same argu- 
ment as above. The QCD corrections for the rest = 0.55, r)z = 0.57 and we 
approximate r}f T ~ r]^ T . 

Since | < |M^|, the mass difference in the B system is 

\5m B \ = 2\M&\, (41) 

and is useful to constrain \ db and the new physics parameters, X db or X db depending 
on the model considered. Long-distance effects are negligible in the B system, and the 
SM calculation yields |<$771b| = 0.49 ±0.16 ps" 1 , to be compared with the experimental 
value 0.489 ±0.008 ps" 1 . 

A second restriction regarding B oscillations comes from the time-dependent asym- 
metry in the decay B° — > tpKs (see for instance Ref. |H21 for a precise definition). This 
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process is mediated by the quark-level transition b — *> ccs and takes place at tree-level, 
with small penguin corrections. The amplitude for the decay can then be written to 
a good approximation as A = AV* b V cs , with A real. Therefore the asymmetry is |H| 
a^pK s = lm\ 1 p Ks , with 

(Mf 2 )* A Mf 2 (Mf 2 )* V cb V* Mg 
^ Ks - |Mf 2 | A \M«\ \Mf 2 \ V* b V cs \M«\ ^ 

and provides a constraint on the combination of phases of B, K mixing and the de- 
cay b — > ccs, which are functions of the CKM CP violating phases and angles. Our 
calculation within the SM gives a^K s = 0.71, and with other choices of parame- 
ters for the CKM matrix the prediction may change in ±0.08. The world average 
is a^K s — 0.734 ± 0.054 [64J. This asymmetry can also be expressed as 

a^ Ks = sin(2/5 + 29 B - 29 K ) , (43) 

with (3 one of the angles of the well-known db unitarity triangle, 



(3 = arg 



VuV t * b 



(44) 



and 9b, 0k parameterising the deviation of the mixing amplitude phases with respect 
to the SM, 



M B M K 

2 ^ = arg 7M^' 2dK = aTg m^- (45) 

In the absence of new physics, or if the extra phases 9 B , &k cancel, a^Ks = sin 2/3. 

The phase of Mf 2 is relatively fixed by the determination of e and Smx- Despite the 
good experimental precision of both measurements, the former has a large theoretical 
uncertainty from and the latter from long-distance contributions. This allows 
9k to be different from zero, but it must be small anyway. The agreement of a^x s 
with experiment then constrains the phase 9b- The asymmetry in semileptonic decays 
depends also on 9b jUS] and does not provide any extra constraint on the parameters 
of these models. 



5.1.4 B® oscillations 



The analysis of B® oscillations is very similar to the previous one for the B° system, 
with 

G%My B B Bs m B2 



127T 2 



(\ t sb y^s,(x t ) + /\ l 



(46) 
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and m B o = 5.370 GeV, f Bs = 230 ± 35 MeV, B Bs = 1.30 ± 0.18 63J. The S (x c ) and 
S (x c , x t ) terms have again been neglected, and the QCD correction for the S (x t ) term 
is rjf t a = 0.55. The new contributions are 

Ak = (\l) 2 VTTSo(xT)+2\i b \lv%S (x t ,x T ), 

A|. = -SX sb Xl bV ^Y (x t ) + ^ V fXl. (47) 

a 

with the terms involving Sq(x c ,xt) in and Yq(x c ) in A Bs discarded. The QCD 
correction factors are the same as for Mf 2 - The SM estimate for |5ms s | = 2|M]| S | is 
17.6 ± 5.9 ps _1 . Experimentally only a lower bound for l^m^J exists, l^m^J > 13.1 
ps _1 with a 95% confidence level (CL), which can be saturated in Models I and I and 
thus provides a constraint not always considered in the literature. Larger values than 
in the SM are also possible. 



5.1.5 D° oscillations 



In contrast with the K° and B° systems, D° mixing is mediated by box diagrams with 
Q = —1/3 internal quarks. This circumstance leads to a very small mass difference 
in the SM, as a consequence of the GIM mechanism. In addition, the approximation 
of vanishing external masses is inconsistent, and with a careful analysis including the 
charm mass an extra suppression ~ (m s /m c ) 2 is found |6lH 16 7j. resulting in |5?Bd| ~ 
10~ 17 GeV. NLO contributions, for example dipenguin diagrams [68J, are of the same 
order, but to our knowledge a full NLO calculation is not available yet. Long-distance 
contributions are estimated to be \Srri£)\ ~ 10~ 16 GeV jSSj . 

On the other hand, the present experimental limit, |<5mc>| < 0.07 ps _1 = 4.6 x 10~ 14 
GeV with a 95 CL, is still orders of magnitude above SM expectations. This limit can 
be saturated in Model I with a tree-level FCN coupling X cu [70 . In Model I with a 
new quark B the GIM suppression is partially removed but we have checked that D° 
mixing does not provide any additional constraint for a mass < 1 TeV. Hence here 
we only discuss Model I. The element M® 2 is then 



M 



G 2 F M^f 2 D B D m D o 

127T 2 



(4* 



where m D o = 1.865 GeV and f D = 215 ± 15 MeV jnj. We assume B D = 1.0 ± 0.3. 
We have omitted the SM terms, whose explicit expression can be found for instance 
in Ref. |HZ|, since they are generically much smaller than the one written above. In 



20 



contrast with K° and B° oscillations, the terms linear in X cu are both negligible due 
to the small masses m s , m^, and we have dropped them 4 . The QCD correction is 
r)z = 0.59. The mass difference is given by |<5m£>| = 2 |M^|, and provides the most 
stringent limit on X cu . 



5.2 K decays 
5.2.1 K+^it+vu 

The importance of the rare kaon decay K + — > 7i + uu in setting limits on the FCN 
coupling X sd has been pointed out before jlHj. This is a theoretically very clean process 
after NLO corrections reduce the scale dependence. The uncertainty in the hadronic 
matrix element can be avoided relating this process to the leading decay K + —>■ TT°e + u 
using isospin symmetry, and then using the measured rate for the latter: 



Bt(K + — > ti + vv) r K +a 



2 



£ \KA L + K dV ?X (x t ) + A K+ 2 . (49) 



2 

l=e,fx,T 



Br(K+ -> 7T°e+u) 2n 2 s^\V us 

The factor = 0.901 [73] accounts for isospin breaking corrections. The charm 
contributions at NLO are [7J X e / L = (10.6 ±1.5) x 10" 4 , X T NL = (7.1 ±1.4) x 10" 4 , and 
the QCD correction to the top term is rjf = 0.994 [7J2]. The function X = Co — 4B 
can be found in Appendix El The top and charm terms have similar size because 
X (x t ) > X l NL but \l d < \ c sd . With Br(K+ -> n°e + v) = 0.0487 we obtain the SM 
value Br(K + — > tt + vv) = (6.4 ±0.6) x 10~ n , where in the uncertainty we only include 
that derived from mt and X l NL , and not from CKM mixing angles. Experimentally 
there are only two K + — > n + uu events [ZO]. The corresponding 90% CL interval for 
the branching ratio is [3.2,48] x 10~ n . 

The new physics contributions are denoted by A^+, and in Models I and II they 
read 

= Cu2zX sd , (50) 

where the factor Cviz in is [j_7J 



7T 2 TVS 



C V2Z = - - ^ = -■ (51) 



V2G F M* 



w 



a 



4 Extending the discussion in Ref. [S3] to the case of D° — D° mixing, we can argue that the 
functions multiplying the terms linear in X cu are the Inami-Lim functions appearing in D° — > 
which are also in this case — Yq |72j . 
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In Model I there is another consequence of the mixing of the top quark not considered 
in these expressions: X tt and X?t are different from unity (hence the function C$ 
corresponding to Z penguins changes) and there are extra penguin diagrams with T 
and t, proportional to the FCN coupling XtT- This is the same kind of modification 
that we have seen in the discussion of the radiative correction to Rb. There it was 
found that the net effect of the top mixing would cancel for ttit = m t and is small for 
irtT ~ m t . The magnitude of the correction required to take this effect into account 
in Br(i^ + — > 7i + uu) can be estimated in analogy with that case. We find that the 
correction grows with X t T and my; however, these cannot be both large, as required 
by oblique parameters. The result is that the error made using Eqs. ()49I50|) is smaller 
than the combined uncertainty from X^l and m t (10%). For X%r in its upper limit 
it amounts to a 6% extra systematic error, unimportant with present experimental 
precision. For each value of X t T and ttit we include the estimate of the correction 
required in the total theoretical uncertainty. Bearing in mind the approximation done 
in using Eqs. ()49l5Uj) . we also omit the QCD factors in the calculation because they 
represent a smaller effect. 

This decay sets relevant constraints on A* d , XJ d and X sd that cannot be obtained 
from the rest of processes studied in this Section. This fact has been explicitly proved 
studying what the range of predictions for By(K + — > ir + vv) would be if the rest of the 
restrictions were fulfilled but not the one regarding K + — > 7i + uu. Since in some regions 
of the parameter space Br(i^ + — > 7i + uu) would be out of the experimental interval, 
this process cannot be discarded in the analysis. 

5.2.2 K L -> 

A complementary limit on A* d , and X sd comes from the short-distance contribu- 
tion to the decay Kl — > Although theoretically this is a clean calculation, 
the extraction from actual experimental data is very difficult. The branching ra- 
tio ~By(K l — > /i + /U~) = 7.18 ± 0.17 x 10~ 9 [THj, can be decomposed in a dispersive 
part [Re A] 2 and an absorptive part [ImA] 2 . The imaginary part can be calculated 
from ~Bt(Kl — > 77) and amounts to (7.07 ± 0.18) x 10~ 9 [Sj, which almost saturates 
the total rate. The extraction of the long-distance component from the real part 
[Re A] 2 = (1.1 ±2.4) x 10~ 10 is not model-independent [7§], but as long as our aim is to 
place limits on new physics we can use the model in Ref. [SO] as an estimate, obtaining 
the 90% CL bound ReA SD < 1.9 x 10" 9 . 
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On the theoretical side, the calculation of the short-distance part of the decay is 
done relating it to K + — > /U + z/, 

Bi(K L -» ) SD _ t Kl a 2 

x V^Re + vlY (x t ) Re A' sd + A Kl ] 2 , (52) 

with tk l = 5.17 x 10~ 8 s, tk+ = 1.238 x 10~ 8 s. The factor in the charm term 
is Y NL = (2.94 ± 0.28) x 1(T 4 at NLO [74 . The function Y = C - B can be 
found in Appendix [Bl and the QCD correction for the top is very close to unity, 
7]J = 1.012 [7S]. Using Br(i^ + — > /i + ^) = 0.6343 from experiment, the SM prediction 
is By(Kl — > /i + yU~)sD = (6.6 ± 0.6) x 10~ 10 . The new physics contributions are 

= ReC U2Z X sd . (53) 

As in K + — ► 7r + z/z/ the mixing of the top quark modifies the Inami-Lim function Co 
and adds a new t — T term. The net effect is small and has been taken into account in 
the theoretical uncertainty. 



5.3 B decays 
5.3.1 B -> A s7 

The inclusive decay width T(B — > X^) can be well approximated by the parton- level 
width F(b — * 57). In order to reduce uncertainties, it is customary to calculate instead 
the ratio 

#7 = 777 54 

and derive T(6 — > sj) from i? 7 and the experimental measurement of T(b — > cez?). The 
ratio _R 7 is given by 

I A* I 2 6<t 

^ = W^7W |C ^" 2 ' (55) 

where 2; = m c /mb (pole masses) and 

f(z) = 1 - 8z 2 + 8/ - z 8 - 24/ log 2 (56) 

is a phase space factor for b — > cev. The Wilson coefficient CV 7 (/i) is obtained from the 
relevant coefficients at the scale Mw by RG evolution. 
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The study of b — > 57 in the context of SM extensions with up and down quark 
singlets has been carried out at leading order (one loop) in Refs. (HU E2] • The Wilson 
coefficients at the scale Mw relevant for this process are, using the notation of Ref. [57j . 



C 2 {M W ) 



Kb 
Kb 



t > 



C 7l (M w ) = -^D' (x t ) + AC 7l {M w ) , 

C 8G (M W ) = ~E' (x t ) + AC 8G (M W ) , 
C 3 (M W ) = AC 3 (M W ), 

C 7 (M W ) = AC 7 (M W ), 

C 9 (M W ) = AC 9 (M W ). (57) 
The extra terms in Model I are straightforward to include, 

AC l 7l {M w ) = -l^D' (x T ), 

Z A sb 

AC\ G (M W ) = -l^E' (x T ), 

1 A sb 

AC\{M W ) = 0, 
AC\{M W ) = 0, 

AC\{M W ) = 0. (58) 

In Model E there are contributions from Z penguins with one or two FCN couplings, 
plus H penguins and other terms originated by the non-unitarity of V. The expressions 
read [HI] 

ACI(M W ) = ^(| + C Z + ef) + ^f^KlM+^B)], 

ac&(m*o = ^(I-3ef-3ef)-3^pt[ei(y B )+efw] 

1 X sb 



ACl{M w ) 



6 A 



t ' 



ACj(M w ) = ~^s w - t 



2 2 X sb 



.sb 



^Cl{M w ) = -(l-s 2 w )-£, (59) 



with i/i = (rrii/Mz) 2 , Wi = (mi/M H ) 2 . The functions £ are given in Appendix |Bj 
We have made the approximations y s — 0, y\, — 0, m s /mb = 0, and a very small 
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term proportional to X^X^b has been omitted in ACy 7 and ACg G . Note also that 
C^Myy) = 1 in the SM, but not necessarily in Models I and E. The RG evolution to 
a scale /i = 5 GeV gives [HH EHj 

C 77 (/i) = -0.158 C 2 (M W ) + 0.695 C 7j (M w ) + 0.085 C SG (M W ) 

+0.143 C 3 (M W ) + 0.101 C 7 (M W ) - 0.036 C 9 (M W ) . (60) 

From this coefficient we get _R 7 = 2.62 x 10~ 3 in the SM. In order to incorporate 
NLO corrections we normalise our LO calculation to the NLO value jHl] with an ad 
hoc factor 5 K 1 = 1.12, and we keep the normalising factor for the calculation of 
i? 7 in Models I and E. This is adequate provided the nonstandard contributions are 
small. The systematic error of this approximation is estimated to be smaller than 
~ K^^Kry — 1)[-R 7 — i? 7 (SM)], and vanishes if the new physics terms scale with the 
same factor fT 7 . We have found that in practice this error is of order O(10~ 5 ), and in 
the worst case 2 x 10~ 4 , smaller than the uncertainties present in the LO and NLO 
calculations. With this procedure and using Br(6 — > cev) = 0.102 we obtain Br(6 — > 
sj) = 3.34 x 10~ 4 , in very good agreement with the world average (3.3 ± 0.4) x 10~ 4 
(EEl IHEl EI| • We take as theoretical uncertainty the one quoted in Ref. [Hj, 0.33 x 10~ 4 . 



5.3.2 B X s l + l 



+ /- 



The analysis of the decay B — > X s l + l~ is very similar to the previous one of B — > X s 7. 
Again, the process can be approximated by b — ► sl + l~ and the quantity theoretically 
obtained is the differential ratio 

Ril(s) = v(h 1 ^ J(b - sl + r) , (61) 
1 (o — > cev) as 

with s = (pi+ +pi~) 2 /ml the normalised invariant mass of the lepton pair. The calcula- 
tion at LO [HE] involves two more operators at the scale Mw, Qw and Qioa- Defining 
for convenience Cgy, Cioa by 

Ol ~ OL ~ 

Cw — ^—Cw, Cioa = — Cioa , (62) 

the latter are 

C 9V (M W ) = -^ + ^A_ AZo (x t ) + AC 9V (M w ), 

y a s6 s w 

C 10A {M W ) = J^pl + AC WA (M W ) . (63) 
$w 

5 We obtain the factor K 1 comparing our LO and the NLO calculation of i? 7 in Ref. [HI] with a 
common set of parameters. 
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The extra terms in Model I are analogous to the ones corresponding to the top, 



AC l 9V (M w ) 
AC{ 0A (M W ) 
and in Model E we have jEHj 

AC*y(M W ) 

AC* 0A (M W ) 



A sb 


Y (x T ) 


A sb 


s 2 


A s 


hY (x T ) 


X* c 2 
A sb b W 



(64) 



Cu2Z X s b 

~~ e 2 W 

b W A sb 



X s b 
U2Z^T- > 
A sb 



(65) 



The RG evolution to a scale fi = 5 GeV gives the coefficients of the relevant operators, 

CiO) = -0.221 C 2 (M W ), 
C 2 (ii) = 1.093 C 2 (M W ), 
C w {ji) = C 9V (M W ) + 1.838 C 2 (M W ), 

C 10A (fj) = C 10A {M W ), (66) 
and CV 7 as in the process b — > 57. We define for brevity in the notation an "effective" 

Cfy(fx) = C w (fx)+g(z,s) [3d(/i) + C 2 (/i)] , (67) 
where the s-dependent function g is 



8, 8 16 

■ - log z H 1 

9 & 27 9 s 



2 
9\ 



4z 2 



4z 2 



2 arctan 



x < 



log 



l-A/l 



+ 7ri 



s < 4z 2 



s > 4z 2 



(68) 



Then, the differential ratio Ru is written as 



Ms) 



a 



l \t 12 



47r 2 /(^) |K b | 2 



1-sl 



;i + 2s)(|c-| 2 + |c 



10A| 



+4(1 + 2/s)|CV 7 | 2 + 12ReC; 7 (7 9 f 



(69) 



The partial width Br(6 — > s/ + / ) is derived integrating s from 4m 2 /m 2 to one and 
multiplying by the experimental value of Br (b — > ceu). Within the SM we obtain 
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Br(6 -> se + e-) = 7.3 x 10~ 6 , Br(b -> s^+^u^) = 5.0 x 10" 6 . These values are a little 
sensitive to the precise value of used. We use, as throughout the paper, the MS 
definition. Experimentally, Br(6 — > s[i + fi~) = (8.9 ± 2.7) x 1CT 6 but for electrons only 
an upper bound exists, Br (b -> se+e~) < 11.0 x 10~ 6 with a 90% CL [H]. The NLO 
values IHIj are 10% larger but in view of the experimental errors it is not necessary to 
incorporate NLO corrections to set limits on new physics. The theoretical uncertainties, 
including the possible modification of the Cq functions, do not have much importance 
either (in contrast with the decay b — > 57) and we do not take them into account in 
the statistical analysis. 

Despite the worse experimental precision, in Model E b — > sl + l~ sets a stronger 
limit on the FCN coupling X s b than b — > 57. This is understood because in this model 
the decay b — > sl + l~ can be mediated by tree-level diagrams involving X s b, while in 
the process b — > 37 this vertex appears only in extra penguin diagrams and unitarity 
corrections, of the same size as the SM contributions. Both I = e and I = // have to 
be considered, as the former sets the best upper bound and the latter provides a lower 
bound. In Model I it also gives a more restrictive constraint than b — > 57, but we still 
include the latter in the fit. 



5.4 The parameter e' 

This parameter measures direct CP violation in the kaon system (its definition can be 
found for instance in Ref. |H2])- F° r several years the experimental measurements have 
been inconclusive, but now the determination has settled, with a present accuracy of 
~ 10 %. On the contrary, the theoretical prediction is subject to large uncertainties. 
Instead of calculating e' directly, we calculate e'/e ~ Ree'/e, using the simplified 
expression jH2j 

^ = F £ ,(x t )ImA<, + A £ ,, (70) 

with 

F £ ,(x t ) = P + P x X (x t )+P Y Y (x t )+P z Z (x t )+P E E (x t ) (71) 

and A e i representing the new physics contribution. The factors multiplying the Inami- 
Lim functions are 

P = -3.167+ 12.409 B { 6 1/2) + 1.262 Bf /2) , 
Px = 0.540 + 0.023 Bi 1/2) , 
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Py 



Pe 



Pz 




(72) 



with Bq ' ', Bf Z) non-perturbative parameters specified below. The new contributions 
are 



In A l £ , we approximate the P coefficients in F £ >(xt) with the corresponding ones in 
F e ,(x t ). From lattice or large N c calculations B^ 1/2) = 1.0 ± 0.3, B*£ ,2) = 0.8 ± 0.2. 
Corrections accounting for final state interactions [HH1 EH] modify these figures to 
B£ /2) = 1.55 ± 0.5, B^ /2) = 0.7 ± 0.2 yielding e'/e = (1.64 ± 0.70) x 10" 3 in good 
agreement with the world average (1.72 ±0.18) x 10~ 3 jHIHj- With large N c expansions 
at NLO very similar results are obtained [US]. Notice that the large theoretical error 
in the B parameters partially takes into account the different values from different 
schemes. These uncertainties, together with cancellations among terms, bring about a 
large uncertainty in the prediction. In spite of this fact e'/e is very useful to constrain 
the imaginary parts of A* d , X^ d and X sd [77]. 

5.5 Summary 

The combined effect of the low energy constraints from K and B physics is to disallow 
large cancellations and "fine tuning" of parameters to some extent. As emphasised 
at the beginning of this Section, the various observables depend on the CKM angles 
Vtd, Vt s , V t b and the new physics parameters in different functional forms. Therefore, 
if theoretical and experimental precision were far better the parameter space would 
be constrained to a narrow window around the SM values, and perhaps other possible 
regions allowed by cancellations. However, as can be seen in Tabled present theoretical 
and experimental precision allow for relatively large contributions of the new physics 
in Models I and E, with large deviations in some observables. 

Other potential restrictions on these models have been explored: the rare decays 
B — ► l + l~, B s — ► l + l~ and B — > X s vv. With present experimental precision they do 
not provide additional constraints, nor the predictions for their rates differ substantially 
from SM expectations. An important exception is By(Kl — > ir°uu). This decay mode 




F £l (x T ) ImAj, 

C U2Z (Px + Py + Pz)lmX sd 



(73) 
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SM prediction Exp. value 



e 


2.18 x 10 


-3 


(2.282 ±0.017) x 10~ 3 


\5m B \ 


0.49 




0.489 ± 0.008 




0.71 




0.734 ±0.054 


\8m Ba \ 


17.6 




> 13.1 (95%) 


\5m D \ 


~ 10" 5 




< 0.07 (95%) 




6.4 x 10" 


11 


[3.2 - 48] x 10- 11 (90%) 


Bi(K L -> /U+//")sd 


6.6 x 10 _ 


10 


< 1.9 x 10- 9 (90%) 


Br(6 — > 57) 


3.34 x 10 


-4 


(3.3 ±0.4) x 10- 4 


Br(6 — > se + e~) 


7.3 x 10" 


-6 


< 11.0 x 10- 6 (90%) 


Br(6 — > s{i + fi~) 


5.0 x 10" 


-6 


(8.9 ±2.7) x 10- 6 




1.64 x 10 


-3 


(1.72 ±0.18) x 10- 3 



Table 5: Experimental values of the low energy observables used in the fits, together 
with the SM calculations with the parameters in Appendix El The theoretical errors 
can be found in the text. The mass differences are in ps _1 . 



does not provide a constraint yet, but in Models I and E it can have a branching ratio 
much larger than in the SM. Its analysis is postponed to Sectional 



6 Other constraints 

The diagonal couplings of the u, d quarks to the Z boson are extracted from neutrino- 
nucleon scattering processes, atomic parity violation and the SLAC polarised electron 
experiment (see Ref. jHj and references therein for a more extensive discussion). The 
values of R , c d L R derived from vN neutral processes have a large non-Gaussian cor- 
relation, and for the fit it is convenient to use instead 
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= arctan 
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= arctan 
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C R 





(74) 
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The SM predictions for these parameters, including radiative corrections [5] and using 
the MS definition of s^-, are collected in Table |Hl together with their experimental 
values. The correlation matrix is in Table 





SM 


riYnpfirnprita 1 

J-JAUCI llllv^lllidl 






prediction 


measurement 


Error 


9l 


0.3038 


0.3020 


0.0019 


9n 


0.0300 


0.0315 


0.0016 


0L 


2.4630 


2.50 


0.034 


Or 


5.1765 


4.58 


+0.40 
-0.27 



Table 6: SM calculation of g\, g R , 9l, 6 r and experimental values. 



A & h e R 



sL i 

9l 0.32 

6 L -0.39 

e R ~ o ~ o 



-0.39 ~ 

-0.10 ~0 

1 0.27 
0.27 1 



0.32 
1 

-0.10 



Table 7: Correlation matrix for the parameters g L , g R , 9l, 9 r . 



The interactions involved in atomic parity violation and the SLAC polarised electron 
experiment can be parameterised with the effective Lagrangian 

G 

£ = — j= \p\iZli>P1 h z QiYQi + C 2 ie^f^e qa^qi (75) 
V 2 i=u ^ d 

plus a QED contribution. We are not considering mixing of the leptons, and hence the 
coefficients Cu, Cn are at tree level [TH] 



X 4 
2~ ~ 3' 



'-'In — \ n n b W I > 



/ < I Xdd , 2 2 

— — I o r ~S W I , 



C*2m — X uu ( — — + 2s^/ I . 



C 2d = -X dd [- l - + 2s 2 w ) . (76) 
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The parameters C\ u and Cid can be extracted from atomic parity violation measure- 
ments. The combination C 2 = C 2u — C 2 d/2 is obtained in the polarised electron exper- 
iment. In TablelHlwe quote the SM predictions of C\ u , Cid, C2, including the radiative 
corrections to Eqs. (|7fij) . and the experimental values. The correlation matrix is in 
Table El 





SM 
prediction 


Experimental 
measurement 


Error 


Clu 


-0.1886 


-0.209 


0.041 


Cid 


0.3413 


0.358 


0.037 


c 2 


-0.0492 


-0.04 


0.12 



Table 8: Experimental values and SM calculations of the parameters C\ u , Cid i n the 
Lagrangian in Eq. (|75|l and the combination C 2 = C 2u — C 2 d/2. 





Clu 


Cm 


c 2 


Clu 


1 


-0.9996 


-0.78 


Cid 


-0.9996 


1 


0.78 


c 2 


-0.78 


0.78 


1 



Table 9: Correlation matrix for the parameters Ci u , Cid and C 2 = C 2u — C 2 d/2. 

The leading effect of mixing with singlets is a decrease of X uu (in Model I) and 
Xdd (in Model E), which reduces g\ and the modulus of C 2 in both cases, and also the 
modulus of Ciu or Cid, depending on the model considered. The angle 6l grows when 
the up quark mixes with a singlet and decreases when the mixing corresponds to the 
down quark. The right-handed couplings are not affected by mixing with singlets, and 
thus g\ and 9r remain equal to their SM values. However, they have to be included 
in the fit because of the experimental correlation. Another consequence of the mixing 
may possibly be the modification of the small radiative corrections to Eqs. ()74|76j) . We 
neglect this subleading effect and assume that the corrections remain with their SM 
values. 
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7 Some observables from K and B physics with 
large new effects 



There is a large number of observables of interest that will be tested in experiments 
under way and for which our models lead to departures from the SM. Necessarily, our 
study is not complete and we pass over many relevant processes that deserve further 
attention. We discuss for illustration the CP-violating decay K° — > tx^vv and the 
time-dependent CP asymmetry in B° s — > D+D~ . 



7.1 The decay K° -> ^vv 

This decay is closely related to K + — > 7i + uu, but its detection is much more difficult. 
At present it is still unobserved, and the 90% CL limit on this decay mode is Br(K L — > 
tx^vv) < 5.9 x 1(T 7 96 . It is calculated as 



Bi{K L -> it°vu) = r KL ^Br(K + -> vr°e + z/ 



x 



r K + 2ir 2 s^ v \V us \ 2 
vfX {x t ) Im A* d + Im A^+l 2 , (77) 



with the corresponding isospin breaking correction tk l = 0.944. The SM prediction 
for this partial width is (2.3±0.17) x 10 -11 , one third of the value for By(K + — > 7i + uu). 
As only the imaginary parts of the CKM products enter in the expression for Bt(Kl — > 
tt°uP), it is possible to have this decay rate much larger while keeping By(K + — > ix+vv) 
in agreement with experiment. However, there is a model-independent limit 

Br(K+ -> vr+z/z/) ~ v ; 

that holds provided lepton flavour is conserved [HZj- We will find later that in our 
models this bound can be saturated, leading to a increase in Bt(Kl — > tt uu) of an 
order of magnitude, even keeping Br(K + — > 7r + z/z/) at its SM value. 



7.2 The CP asymmetry in 5 S ° D+Uj. 

The "gold plated" decay — > D+D~ is mediated by the quark- level transition b — > 
ccs. This transition is dominated by a single tree-level amplitude A cx V* b V cs and it is 
then free of hadronic uncertainties [98J. The time- dependent CP asymmetry is written 
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aS a D+DJ = lmX D+DJ' With 

\M B «\ V*V ^ ' 

\ 1V1 12 I V cb V cs 

and also 

a 0s+i?r =sin(2C-2^J, (80) 

with 



C = arg 



v tb v t * s _ 

B s 



*1) 



and 9b s parameterising the deviation of the phase of M 12 s with respect to its SM value, 



M12 

(A^12~ JSM 



In the SM ( ~ 0, so the asymmetry a D + D - is predicted to be very small (a D + D - ~ 0.03 
with the parameters for the CKM matrix in Appendix El). Therefore, its measurement 
offers a good opportunity to probe new physics, which may manifest itself if a nonzero 
value is observed pi] . 

Another possible final state given by the same quark-level transition is ip<p. This 
state has a clean experimental signature at hadron colliders, ip — > l + l~ and <fi — > K + K~ , 
providing better chances to measure sin(2C — 28b s ) at Tevatron However, in this 

case both particles ip and have spin 1, and then the orbital angular momentum 
is not fixed. (In B° s — > DfD~ the two D mesons have spin 0, and therefore they 
are produced in a / = CP-even state.) The are produced in an admixture 
of CP-even and CP-odd states, which can be disentangled with an analysis of the 
angular distribution of their decay products / + /~, K + K~ |101j . We will loosely refer 
to the CP asymmetries containing the phase (2( — 29 Bs ) as a D + D ~, understanding 
that this includes the asymmetries corresponding to final states DfD~, ifxp, D* + D*~ , 
etc. Finally, it is worthwhile mentioning that cos(2C — 26 b s ) can also be measured, on 
condition that the width difference between the two mass eigenstates is sizeable [102J. 
This can be done without the need of tagging the initial state (B® or B®) and provides 
an independent measurement of this important phase. 

8 Results 

We are more interested in the departures from SM predictions originated by mixing 
with exotic singlets than in finding the best fit to all experimental data. Therefore 
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we must specify the criteria for what we will consider as agreement of these models 
with experiment, and of course our predictions depend on this choice. We require: 
(?) individual agreement of observables with data, and (n) that the joint x 2 °f the 
observables, divided in subsets, is not much worse than the \ 2 of these subsets in 
the SM. For the first condition, the number of standard deviations allowed in a single 
observable is similar to the departure already present in the SM. The second condition 
consists in requiring that the x 2 of a group of variables is smaller or equal to the SM 
X 2 value increased by a quantity numerically equal to the number of variables in the 
group. In average, this condition means admitting a la deviation for a variable which 
in the SM coincides with the experimental measurement, an extra departure of 0.41<r 
for a variable which is at la within the SM, or 0.24a for a variable which is already at 
2a in the SM. This second condition is in practice much stronger than the first one. 
These criteria are best explained by enumerating them: 

• The moduli of the CKM angles V u d, V us , V u b, V c d, V cs , V c b can be at most 2a away 
from the figures in Table d The sum of the x 2 must be smaller than the SM 
result plus six. 

• The predictions for Rb, R c , A F ' B , A F %, Ab, A c may be up to 3a away from the 
central values in Table El We allow larger departures than in the previous case 
because the SM prediction of A F ^ is almost 3a from the experimental measure- 
ment. We also require that the x 2 (calculated with the correlation matrix in 
Table EJ) is smaller than the SM result plus four (in Model I) or plus three (in 
Model E). The number of variables in this subset that effectively change with the 
mixing are four and three in Models I and E, respectively. 

• The contributions to the oblique parameters S, T, U from new physics have to 
be within 2a of the values in Table |U The sum of the x 2 °f the three variables 
must be smaller than the SM x 2 value plus three. 

• The observables e, \5i71b\ and a^x s are allowed to move within 2.5a of the num- 
bers in Table El The total x 2 has to be smaller than the SM value plus three. 

• The branching fractions for b — > 57 and b —* s[i + are required to agree with 
the experimental figures in Table El within 2a, and their x 2 has to be smaller than 
the SM result plus two. 

• The departure of e'/ e from the experimental measurement in Table El can be at 
most la larger than the departure within the SM. 
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• The parameters g 2 Ll g R , 9l, Or, C\ u , and C 2 have to be within 2a of the central 
values in Tables El and |H1 Their \ 2 computed with the correlation matrices in 
Tables [3 and El is required to be smaller than the SM value plus five (gR and Or 
are not affected by the mixing in these models). 

The observables e, {Smsl, Br(6 — > 57) and e'/e have large theoretical errors that are 
of similar magnitude as the experimental ones. In the comparison of these observables 
with experiment we use the prescription explained in Appendix [0 assuming that the 
theoretical errors are Gaussian. For |5t71d|, Br(K + — ► tt + uu), Bt(Kl — ► /U + /U~)sd arid 
Br (b —>■ se + e~) we require that the predictions are in the experimental intervals quoted 
in Table El (the upper limit of |5m£>| in the literature has a 95% CL, instead of the 
more common one of 90%). We also set the condition \5mB a /SmB\ > 26.7, rather 
than |5?71bJ > 13.1 ps -1 , to avoid theoretical uncertainties. With all these restrictions 
we explore the parameter space to find the interval of variation of charged current 
top couplings, flavour- diagonal and flavour- changing Z couplings and the observables 
introduced in last Section. We discuss the results for Models I and E separately. 

8.1 Mixing with an up singlet 

One fundamental parameter in Model I is the mass of the new quark %. We find 
that for low itlt the effects of mixing can be huge, with Vtd, Vt s and Vtb very different 
from the SM predictions. In this scenario the effects of mixing on Rf, and oblique 
corrections almost cancel, while the new quark can virtually take the place of the top 
in reproducing the experimental values of the meson observables analysed. This is 
viable because for small tjit the Inami-Lim functions for the top and the new quark 
are alike. As itlt grows this possibility disappears and the dominant contributions 
come from the top, but still significant departures from the SM can be found. In our 
analysis we have checked that for the values and plots shown the Yukawa couplings 
remain perturbative. The decoupling limit is not reached in any of the cases considered. 

One of the most striking results in Model I is the deviation of \V t b\ from unity (see 
Fig. H}. The modulus of Vtb is determined by the coupling Vrb in Fig. El and the latter 
is bounded by the T parameter, as it was seen in Section |U (The dependence on only 
one observable leads to the very simple behaviour of the curves in Figs. QancEl) For 
ttlt = 200 GeV, \V t b\ can be as small as 0.58. The lower limit on \V t b\ grows with %, 
but even for m T = 600 GeV it is \V t b\ > 0.977, substantially different from the SM 
prediction \Vtb\ = 0.999. Although sizeable and theoretically very important, this 2% 
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difference is difficult to detect experimentally at LHC, which is expected to measure 
the size of V t b with a precision of ±0.07 |13j . 




Figure 1: Allowed values of \V t b\ (shaded area) in Model I, as a function of the mass of 
the new quark. 
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Figure 2: Allowed values of the coupling \ Vrb\ of the new quark (shaded area) in Model 
I, as a function of its mass. 

The top charged-current couplings Vt s and Vtd can be very different from SM ex- 
pectations as well. In the SM 3x3 CKM unitarity fixes \V ts \ — \V c b\. In Model I \Vf S \ 
can be between 0.002 and 0.061 for m T = 200 GeV (see Fig. EJ). The allowed interval 
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narrows as my increases, and for my = 600 GeV the interval is essentially the same as 
in the SM. The range of variation of V t d is also considerably greater than in the SM 
(see Fig. HJ). For uit < 300 GeV Vtd can be almost zero (and in this case the T quark 
would account for the measured values of K and B observables) , or even larger than 
Vtsi as can be seen in Fig. EJ Again, for heavier T the permitted interval decreases and 
for tjit = 600 GeV it is practically the same interval as in the SM. We remark that 
the curves in Figs. OH giving the upper and lower bounds arise from the various re- 
strictions discussed in Sections OKI especially those regarding meson observables, thus 
their complicated behaviour should not be surprising. We do not claim that the blank 
regions in these three figures are excluded. The quoted allowed limits might be wider 
if some delicate cancellation not found in the numerical analysis allows a small region 
in parameter space with Vtd, V ts or their ratio outside the shaded areas. 




600 



Figure 3: Allowed values of \V ts \ (shaded area) in Model I, as a function of the mass 
of the new quark. 

In contrast with the former, the intervals for CKM mixing angles Vrd, Vr s do not 
show a pronounced decrease with m^. Vrd can be in the interval < \Vrd\ < 0.05 for 
the rriT values studied, and the maximum size of \Vt s \ decreases from 0.06 for my = 200 
to 0.05 for m T = 600 GeV. 

The counterpart of the departure from the SM prediction \V t b\ = 0.999 is the de- 
crease of the Ztt coupling. Within the SM, the isospin-related term X u equals one by 
the GIM mechanism, while in Model I the GIM breaking originated by mixing with a 
singlet reduces its magnitude. The modulus of X tt , as well as V#, is determined by the 
parameter Vrb and hence its possible size is dictated only by the T parameter. The 
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Figure 4: Allowed values of \V t d\ (shaded area) in Model I, as a function of the mass 
of the new quark. 




m T 

Figure 5: Allowed values of the ratio |Vtd/^tsl (shaded area) in Model I, as a function 
of the mass of the new quark. 

interval allowed for X tt is plotted in Fig. where we observe that for = 200 GeV 
it reaches down to X tt = 0.34. The lower limit of the interval grows with mj- and 
is approximately X u = 0.96 for tjit = 600 GeV. The Ztt coupling will be precisely 
measured in ti production at TESLA. With a CM energy of 500 GeV and an integrated 
luminosity of 300 fb _1 , 34800 top pairs are expected to be collected at the detector in 
the semileptonic channel Ivjjjj, with / an electron or a muon. The estimated preci- 
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sion in the determination of X tt with this channel alone is of 0.02. Then, even with 
ttlt = 600 GeV a 2a effect could be visible. 

l 

0.9 
0.8 

_ 0.7 

>< 

0.6 
0.5 
0.4 

200 300 400 500 600 

m T 

Figure 6: Allowed values of \X tt \ (shaded area) in Model I, as a function of the mass 
of the new quark. 

FCN couplings are perhaps the most conspicuous manifestation of mixing with 
quark singlets, and offer another excellent place to search for new physics. In the SM 
they vanish at tree-level by the GIM mechanism, and the effective vertices generated 
at one loop are very small as a consequence of the GIM suppression |103j . This results 
m a neglig ible branching ratio Br(t -> Zc) ~ 10~ 14 within the SM. In Model I the FCN 
coupling X ct can be sizeable , leading to top decays t —> Zc |104j , Zt production at 
LHC dHUHEl and single top production at linear colliders [Ml CHHl IM ED] • For 
irtT ~ m t the new contributions to meson observables involving T diagrams are small, 
and this FCN coupling can be relatively large, \X ct \ = 0.036 (see Fig. Ej) 6 . A coupling 
of this size yields a branching ratio Br(t — > Zc) = 6.0 x 10 -4 (nine orders of magnitude 
above the SM prediction) that would be seen at LHC with 18cr statistical significance 
in top decays and 4.6a in Zt production (with an integrated luminosity of 100 fb _1 ), 
and at TESLA with 8.2a significance in single top production (with 300 fb _1 ). For 
larger itlt, the contributions of the T quark to meson observables (in particular to 
K + — > -n + vv and the short-distance part of — > yU + /i~) decrease monotonically the 
upper limit on \X ct \, with some very small local "enhancements" that can be observed 

6 The reduction with respect to the number quoted in Ref. jHl] is mainly due to the improved limit 
on |<5m£>|. 
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in Fig. [7| For T very heavy there is still the possibility of \X ct \ = 0.009, giving 
Br(£ — > Zc) = 3.8 x 10~ 5 , which would have a 1.2a significance in top decay processes 
at LHC. 




Figure 7: Allowed values of the FCN coupling \X ct \ (shaded area) in Model I, as a 
function of the mass of the new quark. 

In Model I the X ut coupling can have the same size as X ct . This contrasts with other 
SM extensions (for instance, SUSY or two Higgs doublet models) where observable FCN 
tc vertices can be generated but tu vertices are suppressed. The observability of a Ztu 
FCN coupling is the same, and even better in the case of Zt production processes 
at LHC. The coupling X tT between the top and the new mass eigenstate (which is a 
function of the charged-current coupling Vrb) can reach the maximum value permitted 
by the model, X t r = 0.5 for tjit < 210 GeV, descending slowly to a maximum of 
X tT = 0.2 when m T = 600 GeV. 

The mixing with a new singlet may also give new effects in low energy observables. 
The branching ratio of Ki — ► tt°uu can reach 2 x 10~ 10 for "low" my, and 4.4 x 10~ 10 
for mj- > 300 GeV, one order of magnitude above the SM prediction By(Kl — > ifivv) = 
2.4 x 10 -11 . These rates would be visible already at the E391 experiment at KEK, which 
aims at a sensitivity of 3 x 10~ 10 , and up to ~ 10 3 events could be collected at the 
KOPIO experiment approved for construction at BNL (for a summary of the prospects 
on the rare decays K + — > tt + vv and Kl — > ^vv see for instance Ref. |lllj ). The ratio 
Br(i^^ — > 7r z/z/)/Br(i^ + — > tt + vv) of the decay rates of the two kaon "golden modes", 
plotted in Fig.^l can be enhanced an order of magnitude over the SM prediction ~ 0.35, 
and saturate the limit in Eq. (|7%j) for itlt > 310 GeV. This enhancement and a larger 
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Figure 8: Allowed values of the coupling \Xfr\ of the new quark (shaded area) in Model 
I, as a function of its mass. 

value of Br(K + — > -k + vu) (compatible with experimental data) lead to the maximum 
value Bi(K L — > ifivv) = AA x 1CT 10 . On the other hand, a strong suppression of 
this decay mode is possible, with values several orders of magnitude below the SM 
prediction. 




Figure 9: Range of variation of By(Kl — > 7r°z/z/)/Br(i ; r + — > it + vv) in Model I (shaded 
area) . 

The mass difference in the B s system is predicted to be ~ 18 ps" 1 within the SM. 
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The existing lower bound |5m#J > 13.1 ps _1 can be saturated in practically all the 
interval of tut studied. The ratio \SmB a /SmB d \ has been proposed for a determination 
of |Vts/Vtd| jS]. Of course, this determination is strongly model-dependent, because 
new physics may contribute to both mass differences. This ratio equals 36 in the SM, 
and in Model I it may have values between the experimental lower limit of 26.7 and 
77. Finally, the asymmetry a D + D -, which practically vanishes in the SM, provides a 
crucial test of the phase structure of the CKM matrix. The non-unitarity of the 3x3 
CKM submatrix and the presence of extra CP violating phases in Model I allow the 
asymmetry a D + D - to vary between —0.4 and 0.4 for % > 275 GeV, as can be observed 
in Fig. EU 

l 

0.5 



V 



-0.5 
-1 

200 300 400 500 600 

m T 

Figure 10: Range of variation of a D + D - in Model I (shaded area). 
8.2 Mixing with a down singlet 

In Model I the mass of the new quark does not play an important role in the constraints 
on the parameters of the model. The only dependence on tub appears in the D° 
mass difference (which at present does not imply any restriction at least for masses 
up to 1 TeV), b — > sj (less restrictive than b — > sl + l~) and oblique parameters, 
which are less important than and have no influence in practice. Agreement of 
the latter with experiment requires that |Vt&| is very close to unity, \V t b\ > 0.998. This 
is indistinguishable from the SM prediction |Vt&| = 0.999, and forces V t d and V ts to 
be within the SM range, 0.0059 < \V td \ < 0.013, 0.035 < \V ts \ < 0.044. The CKM 
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matrix elements involving the new quark are all small, \V u b\ < 0.087, \V c b\ < 0.035, 
WtB\ < 0.041, but noticeably they can be larger than V u b- 

FCN couplings between the light quarks are small (as required by low energy ob- 
servables), especially the coupling between the d and s quarks, \X ds \ < 1.0 x 10~ 5 . It 
makes sense to study Re X ds and Im X ds separately, even though in principle X ds is not 
a rephasing-invariant quantity. This is so because Eq. ()38j) assumes a CKM parame- 
terisation with V* d V us real. This requirement eliminates the freedom to rephase X ds 
(up to a minus sign) and enables to separate its real and imaginary parts meaningfully. 
The region of allowed values for X ds is plotted in Fig. ^2 for comparison with other 
analyses in the literature jHHl EH EE] • This figure must be interpreted with care: the 
density of points is not associated to any meaning of "probability", but it is simply 
an effect related to the random generation and CKM parameterisation used to obtain 
the data points, and the finiteness of the sample. The height of the allowed area is 
determined by the e'/e constraint, and the width by Kl — > Comparing this plot 

with the ones in Refs. |3T)| l3*%] we see that the left part of the rectangle determined by 
e'/e and Kl — > /U + ' n~ is practically eliminated by the constraints from K + — ► n + uu 
and e, except the upper left corner. The height of the rectangle is also smaller, meaning 
that in our case the requirement from e' j e (using the prescription in Appendix |HJ) is 
more stringent. 
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Figure 11: Allowed region for the real and imaginary parts of the FCN coupling X ds 
in Model II . 

The upper bounds for X d b and X s b found in our analysis are \X d b\ < 1.1 x 10~ 3 , 
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\X ab \ < 1.1 x 10~ 3 . Plots analo gous to Fig. ^2 are not meaningful for these parameters, 
because there is a freedom to rephase the b field and change arbitrarily the phases of 
Xdb and X s f). The only meaningful bounds are hence the limits on their moduli. The 
FCN coupling X^b is not so limited by low energy measurements, and can reach 0.041. 

Despite these restrictions on X^s, X^b, X s b and the fact that CKM matrix elements 
involving the known quarks must be within the SM range, the presence of tree-level 
FCN couplings has sizeable effects on some low energy observables, of the same mag- 
nitude as in Model I. A decay rate Bt(Kl — > tt^vv) = 1.6 x 10~ 10 can be achieved, and 
the ratio By(Kl —>■ tt° v v) / 'Br (K + — > 7i + uu) can equal 4.34. The lower limit on |5msJ 
can be saturated, and the ratio \5rriB s /5mB d \ can be up to 67. On the other hand, 
the asymmetry a D + D - takes values between —3 x 10~ 3 and 0.11. This upper limit is a 
factor of 2 — 3 larger than the SM prediction. 



8.3 Mixing with two singlets 



Mixing with more than one singlet lets two quarks of the same charge, for instance the 
d and s quarks, mix significantly with exotic quarks without necessarily generating a 
FCN coupling X^ s between them, in virtue of Eqs. El This allows a better fit to the 
measured CKM matrix elements and u, d diagonal couplings to the Z boson, especially 
in Model E. In this model the global fit can be considerably better than in the SM, for 
instance with the CKM matrix 
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(83) 



The actual masses of the two extra down quarks are not very relevant, and have been 
taken as 200 and 400 GeV in the calculation. In this example the x 2 of the six measured 
CKM matrix elements is 1.14, while in the SM best fit it is 4.77. The parameters 
describing the Zuu and Zdd couplings are g\ = 0.3024, 9i = 2.4612, C\d = 0.3398, 
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C*2 = —.0492, and the others unchanged with respect to their SM values. The x 2 °f 
these parameters is 7.73, improving the SM value of 10.5. The agreement of the rest 
of the observables with experiment is equal or better than within the SM, as can be 
seen in Table El (the experimental results for e and 8m B can be accommodated with 
slightly larger B parameters). This example of "best fit" matrix gives the predictions 
\X bB \ = 0.006, \X hB ,\ = 0.004. The result Br{K L -> tc uu) = 2.0 x 10~ 10 is very 
similar to the SM case, but other examples with a little worse \ 2 can be found, having 
enhancements (or suppressions) of this rate by factors up to three. These examples 
show explicitly that new physics effects are not in contradiction with good agreement 
with experimental data, although our restrictive criteria for agreement with experiment 
at the beginning of this Section already made it apparent. 



Observable 


Value 


Rb 


0.21590 


R c 


0.1724 


A o,b 

^FB 


0.1039 


A 0,c 
^FB 


0.0744 


A 


0.935 


Ac 


0.669 


£ 


2.08 x 10- 


\8m B \ 


0.45 


\5m Bs \ 


17.6 




0.74 



Br(AT+ 


— > 71 + UU) 


6.0 x 10" 11 


Br(K L - 


+ / i+ A i ")sD 


6.3 x lO" 10 


Br(6 


-> si) 


3.35 x 10- 3 


Br(6 - 


■> se + e~) 


7.3 x 10" 6 


Br(6 - 




5.0 x 10~ 6 


£ 


'/e 


1.6 x 10-3 



Table 10: Values of some observables for the "best fit" matrix in Model II with two 
extra singlets. The mass differences are in ps _1 . 

Finally, we have also noticed that the predictions for the parameters and observables 
under study do not change appreciably neither in Model I nor in Model I when we 
allow mixing with more than one singlet of the same charge. 
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9 Conclusions 



The aim of this paper has been to investigate how the existence of a new quark sin- 
glet may change many predictions of the SM while keeping agreement with present 
experimental data. In Model I the mixing with a Q = 2/3 singlet might lead to huge 
departures from the SM expectation for the CKM matrix elements Vtd, V ts , V t b and 
the diagonal coupling Ztt. Additionally, observable FCN couplings Ztu and Ztc may 
appear. These effects depend on the mass of the new quark, as has been shown in 
Figs. Qini For itit ~ 200 GeV the new quark might effectively replace the top in 
reproducing the experimental observables in K and B physics, allowing for values of 
Vtd and Vts very different from the SM predictions. On the other hand, for larger rnx 
the leading contributions to K and B observables are the SM ones, with possible new 
contributions from the new quark. This effect can be clearly appreciated in Figs. 0H3 
where it is also apparent how important a direct determination of Vtd and Vt. s would be. 
Unfortunately, the difficulty in tagging light quark jets at Tevatron and LHC makes 
these measurements very hard, if not impossible. Any experimental progress in this 
direction would be most welcome. 

The mixing of the top with the new quark results in values of Vtb and the Ztt 
coupling parameter X u significantly smaller than one. These deviations from unity 
would be observable at LHC [T3] and TESLA, respectively. For larger m T , \V t b\ and 
\X a \ must be closer to unity, as can be seen in Figs. ^ and H3 However, the decrease in 
X tt would be visible at TESLA even for rriT = 600 GeV. The FCN couplings Ztu and 
Ztc could also be observed at LHC for a wide range of % |104t 1105} 1106] . 

The effects of top mixing are not limited to large colliders. Indeed, the observables in 
K and B physics studied here provide an example where these effects do not disappear 
when the mass of the new quark is large. We have shown that the predictions for the 
decay Kl — > 7r°z/z/, the 5mg s mass difference and the CP asymmetry a D + D - can be 
very different from the SM expectations, and effects of new physics could be observed 
in experiments under way or planned. These predictions for Model I are collected in 
Table ITT1 Before LHC operation, indirect evidences of new physics could appear in the 
measurement of CP asymmetries at B factories. A good candidate is the asymmetry 
a D + D - discussed here, but many other observables and CP asymmetries are worth 
analysing. If no new physics is observed, further constraints could be placed on CP 
violating phases. 

In Model E the effects of the new Q = —1/3 singlet on CKM matrix elements are 
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Quantity 



Range 



\Vtb\ 


0.58 


1 


\Vu\ 


4 x 1(T 5 


0.044 


\Vts\ 


0.002 


0.06 


\v td /v ts \ 


6 x lO" 4 


2.9 


\v Td \ 





0.052 


\Vts\ 





0.063 


\V Tb \ 





0.81 


\X tt \ 


0.34 


1 


\x ut \ 





0.038 


\Xct\ 





0.036 


\X t r\ 





0.5 


Bl(K L 7T°Z/Z/) 

Bt(K l ^tt°uP) 
Br(K+^n+vu) 

\5m B J5m B \ 


~ 
~ 
26.7 


4.4 x 10 
4.35 
77 


a D+D7 


-0.4 


0.4 



-10 



Table 11: Summary of the predictions for Model I. 



negligible and FCN couplings between known quarks are very constrained by experi- 
mental data. However, the predictions for meson observables, summarised in Table \12\ 
are rather alike. In addition, we have shown how the mixing with two singlets can 
improve the agreement with the experimental determination of CKM matrix elements 
and Zuu, Zdd couplings. This can be done keeping similar and in some cases better 
agreement with electroweak precision data and K and B physics observables. 

All the effects of mixing with singlets described are significant, but of course the 
decisive evidence would be the discovery of a new quark, which might happen at LHC 
or even at Tevatron, provided it exists and it is light enough. In this case, the pattern of 
new physics effects would allow to uncover its nature. Conversely, the non-observation 
of a new quark would be very important as well. If no new quark is found at LHC, the 
indirect constraints on CKM matrix elements and nonstandard contributions to meson 
physics would considerably improve. 
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Quantity 


Ranj 


;e 




\Vtb\ 


0.998 


1 




\Vtd\ 


0.0059 


0.013 




\Vts\ 


0.035 


0.044 




\v uB \ 





0.087 




\Vcb\ 





0.035 




\v tB \ 





0.041 




\Xds\ 





1.0 x 10- 


-5 


ReX ds 


-1.0 x 10~ 5 


3.4 x 10" 


-6 


lmX ds 


-2.7 x 10~ 6 


2.4 x 10" 


-6 


\Xdb\ 





1.1 x 10- 


-3 


\X s b\ 





1.1 x 10" 


-3 


\Xbs\ 





0.041 




Br(K L -> ir°vv) 

Br(K L ^n°uD) 
Bt(K+->tt+uu) 

\5m B J5m B \ 


~ 
~ 
26.7 


1.6 x 10" 
4.34 
67 


•10 


a D+D~ 


-3 x 10~ 3 


0.11 





Table 12: Summary of the predictions for Model E. 



A Common input parameters 

Unless otherwise specified, experimental data used throughout the paper are taken 
from Refs. [H E] . We use the results in Ref . |112j to convert the pole masses rrii to the 
MS scheme and to perform the running to the scale Mz- The results are in Table IT3l 
For u, d, s we quote the MS masses at 2 GeV instead of the pole masses. The numbers 
between brackets are not directly used in the calculations. 

The running masses m c (m c ) = 1.28, rfibijfib) — 4.19 are also needed. The lepton pole 
masses are m e = 0.511 MeV, m M = 0.105 and m T = 1.777 GeV. We take M z = 91.1874, 
Y z = 2.4963, M w = 80.398 and M H = 115 GeV. The electromagnetic and strong 
coupling constants at the scale Mz are a = 1/128.878, a s = 0.118. The sine of the 
weak angle in the MS scheme is s\ = 0.23113. 

The CKM matrix used in the context of the SM is obtained by a fit to the six 
measured moduli in Table HJ and is determined by \V US \ = 0.2224, \V U b\ = 0.00362, 
\Vcb\ = 0.0402, and the rest of the elements obtained using 3x3 unitarity. The phase 



48 



rrti 



mArrii 



mi(M z 



m u 
m d 
m c 
m s 
m t 
m b 



(0.003) 
(0.006) 
1.5 
0.12 
174.3 
4.7 



1.22 

164.6 
4.12 



0.0016 
0.0033 

0.68 
0.067 

175.6 

2.9 



Table 13: Quark masses (in GeV) used in the evaluations. The uncertainty in m t is 
taken as ±5.1 GeV. For u, d, s we write the MS masses m,i{2 GeV) instead of the pole 
masses. 



5 in the standard parameterisation [5] is determined performing a fit to e, e'/e, a^K s 
and |<5ms| with the rest of parameters quoted, and the result 5 = 1.014 is very similar 
to the one obtained in the fit in Ref. 0. 



B Inami-Lim functions 



In this Appendix we collect the Inami-Lim functions used in Section |S| The box 
functions F and So appear in meson oscillations. D' , Eq and E' are related to photon 
and gluon penguins. The functions Xq and Yq are gauge-invariant combinations of 
the box function B and the Z penguin function C , X = C — 4B , Y = C — B . 
The function Z is a gauge-invariant combination of photon and Z penguins. Their 
expressions read [HU EZ| 

_, . N 2, Xi(18 - llxi -xj) x?(15 - IQxi + 4x1) , fOA . 

_j_ K™2 f 7rp. Q /v»3 _i_ o T 2 

rf 3 P. />-» 2 O ,y> , '2 rp 2 

^ = - 4(1 -> + 2(1^ 108 <86> 
. . 4 Tx^Xj 4 -|- X{Xj 2 1 

4(1 — Xi)(l — Xj) 4(1 - Xj) 2 (xi - Xj) 

4 — 8Xj + XiXj 2 i /C7\ 

4 (1 -Xj) 2 (xj -Xi) 1 

„ , \ 3a;,- OC n X^X j (x„ I i 

4(xi - - 1) 4{xi - l) 2 {xi - xj) 
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S (xi) 
Z (xi) 



j (^Xj 8^2/ j I 4 
4(xj - l) 2 {xj - Xi) 
Axi — llxf + x 3 



logx 



J ' 



3x? 



4(1 -xO 2 2(1 -x 



■ \3 



log Xi 



(88) 
(89) 



_ 1 18x 4 - 163a:? + 259x 2 - 108x; 

~9 gXl+ 144(x, - l) 3 

32xf - 38x 3 - 15x 2 + I8xi , 
+ ^ tt^ 1 log^, 



5 (x 4 ) 
C (xi) 
X (xj) 
^o(^) 



1 

4 

~8 

X j 



72( Xi - l) 4 



\ X i. 



log Xi 



Xi - 6 3xj + 2 

7 + 7 Tv7 lo S x i 

Xi - 1 



+ 



(z< - I) 2 
3x,- — 6 



,X ; 



X,; -4 



[Xi - l) 2 



X?; 1 [%i 1 



log ^ 
2 log a* 



(90) 
(91) 
(92) 
(93) 
(94) 



The functions appearing in the Z FCNC penguins involved in the calculation of b — > 57 
are [EH 

1 

,Z / o /i „2 



-3 + 2s^) , 



= ^(-3-4^, 



+ 



8 - 30y B + 9j/| - 5j/| 

144(1 -y B ) 3 ' 8(1 -y B ) 4 
16w B - 29w 2 B + 7w% -2w B + Sw% 



+ 



144(1 - w B ) 3 24(1 - w B Y 

where we have approximated y s = 0, y\> = and m s /m;, = 0. 



logt/B , 

logws 



(95) 
(96) 
(97) 

(98) 



C Statistical analysis of observables with theoreti- 
cal uncertainty 

The most common situation when comparing a theoretical prediction xt with an ex- 
perimental measurement x e is that the uncertainty in the former can be ignored. This 
does not happen for some observables analysed in this article, which are subject to low 
energy QCD uncertainties. For example, if we have for e x e = (2.282 ± 0.017) x 10~ 3 
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and x t = (2.42 ±0.42) x 10 -3 , how many standard deviations is x t from x e l To answer 
naively that it is at 8.1a is clearly wrong, and the comparison between both should 
weigh in some way the error on xt- Here we explain how we obtain in such cases a 
reasonable estimate of the agreement between the theoretical and experimental data. 

Let us recall how x e and x t are compared when the former has a Gaussian distri- 
bution with mean fi e and standard deviation a e and x t is error-free and equals fi t (see 
for instance Ref. [US]). The x 2 value is defined as 

X 2 =(^)\ (99) 
and from it the P number is computed as 



P= / f(z;l)dz, (100) 
where f(z; n) is the x 2 distribution function for n degrees of freedom, 

n/2-l g -z/2 

f( z - n) = -— , . . (101) 
J v ' ; 2 n / 2 T(n/2) v ; 

The P value is the probability to obtain experimentally a x 2 equal or worse than the 
actual one, that is, a result equal or less compatible with the theory. Performing the 
integral in Eq. (jlOOj) . 

P^-erf^i-erf^, (102) 

with erf the well-known error function. The probability to obtain an equal or better 
result is 1 — P. For instance, with |/x e — fi t \ = a e we have 1 — P = erf(l/\/2) = 0.68, 
corresponding to one Gaussian standard deviation, as it obviously must be. 

When Xt is not considered as a fixed quantity \xt but has some distribution func- 
tion g(xt) (that may be Gaussian or may not), we use the probability law P(A) = 
J2i P(A\Bi)P(Bi), with J2iP(Bi) — 1> to convolute the x t -dependent P number with 
9- 

P\nt-^x t g(xt) dx t = l- / erf — -= g(x t ) dx t . (103) 

The assumption that Xt = fit without error can be translated to Eq. ()103jl choosing the 
"distribution function" g(x t ) = S(x t — fit), in which case we recover Eq. ()102|) . 

An adequate (but not unique) choice of the function g(x t ) may be a Gaussian. One 
source of systematic uncertainties is often due to the input parameters involved in the 
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theorical calculation (m t , a s , CKM mixing angles, etc.), whose experimental values are 
given by a Gaussian distribution. It is then likely that the distribution function g{x t ) 
(where Xt is also a function of its input parameters) has a maximum at fit and falls 
quickly for increasing \xt — lh\- This feature can be implemented in a simple way by 
choosing g(x t ) as Gaussian, and we expect that the results are not very sensitive to 
the precise shape of the function g(x t ). 

Let us then assume that g(xt) is a Gaussian with mean fit and standard deviation 
<r f . Intuitively, we expect that if a t <^ a e , Eq. ()103j) should reduce to Eq. (jl02J) . This 
is easy to show. Writing the explicit form of g(x t ), 

i x t-IH) 2 

P = 1 - / + °° erf \^~ Xt \ e dx t . (104) 

The limits of this integral can be taken as fit — n at, fit + n &t, with n > 4. The integral 
is negligible out of these limits due to the exponential (the error function takes values 
between and 1). Changing variables to A t = x t — fit, we observe that |A t | <C o e 
under the assumption that a t ^ a e . Expanding the error function in a Taylor series 
to order A t , the integral can be done analytically, 

P = i _ erf " er f ^ . (105) 



y/2 y/2 



For n > 4, erf n/\/2 ~ 1 to an excellent approximation and we obtain Eq. (|1()2|1 . as we 
wanted to prove. 

Results for P values can be expressed in a more intuitive form as standard "number 
of sigma" n a inverting Eq. (|1()2|1 . 



V2err\l- P) , (106) 



with erf -1 the inverse of the error function. However, this n a does not retain the 
geometrical interpretation of the distance between /x e and fit in units of cr e that has 
when o t ~ 0. 

We apply this procedure to the example at the beginning of this Appendix, with 
x e = (2.282 ± 0.017) x 10~ 3 and x t = (2.42 ± 0.42) x 10~ 3 . Assuming for simplicity 
that the distribution of Xt is Gaussian, we obtain the much more reasonable result of 
n a = 2.25. This number must be compared with n a = 8.1, obtained without taking 
into account the theoretical error, i. e. calculating naively {fi t — /i e )/a e . The use of the 
theoretical error in the statistical comparison mitigates the discrepancy and implements 
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numerically, in a simple but effective way, what one would intuitively expect in this 
case. The result n a = 2.25 reflects the fact that the theoretical and experimental values 
can have a good agreement if x% is smaller than its predicted value fit = 2.42 x 1CT 3 
but a bad one if Xt is larger, what is also possible because the theoretical error is of 
either sign. The prescription presented here has also one very gratifying property: if 
we change [i t — > l^t = A** + &Ht, with 6fi t <C cr t , the P value is hardly affected. For 
fi' t = 2.35 x 1CT 3 , n a changes only to 2.23, while the pull calculated naively decreases 
to 4. 

This construction can be generalised when not all the range of variation of x e , x t is 
physically allowed. We write without proof the expression for P in this case. Assuming 
that the physical region is x e > 0, x t > 0, 

r i -i 

//,„ 

~~ L V2a e \ 

max{/x e , \fj, e -x t \} \ 
+ erf -= g(x t ) dx t . (107) 




V2 



Finally, notice that the expressions in Eqs. (jl03H07|) for the P number are not 
symmetric under the interchange of theoretical and experimental data, even if g(fit) is 
Gaussian. This reflects the fact that P(data| theory) ^ P(theory|data), but they are 
related by Bayes' theorem. 
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